arXiV: astro-ph/0112044v 1 3 Dec 2001 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN style file vl.4) 


Hard X-ray states and radio emission in GRS 1915-|-105 

M. Klein-Wolt^*, R. P. Fender^ G. G. Pooley^, T. Belloni^, S. Migliari^’^ 

E. H. Morgan^, M. van der Klis^ 

^ Astronomical Institute ‘Anton Pannekoek’, University of Amsterdam, and Center for High Energy Astrophysics, Kruislaan 403, 
1098 SJ, Amsterdam, The Netherlands 

^ Mullard Radio Astronomy Observatory, Cavendish Laboratory, Madingley Road, Cambridge CBS OHE 
® Osservatorio Astronomico di Brera, Via. E. Bianchi 46, 1-23870, Merate, Italy 
Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 0213, U.S.A. 


1 February 2008 


ABSTRACT 

We have compared simultaneous Ryle Telescope radio and Rossi X-Ray Timing 
Explorer X-ray observations of the galactic microquasar GRS 1915-1-105, using the 
classification of the X-ray behaviour in terms of three states as previously established. 
We find a strong (one-to-one) relation between radio oscillation events and series of 
spectrally hard states in the X-ray light curves, if the hard states are longer than ~ 100s 
and are “well separated” from each other. In all other cases the source either shows a 
low-level or a high level radio emission, but no radio oscillation events. During intervals 
when the source stays in the hard spectral state for periods of days to months, the 
radio behaviour is quite different; during some of these intervals a quasi-continuous 
jet is formed with an almost flat synchrotron spectrum extending to at least the 
near-infrared. Based on the similarities between the oscillation profiles at different 
wavelengths, we suggest a scenario which can explain most of the complex X-ray : radio 
behaviour of GRS 1915-1-105. We compare this behaviour with that of other black hole 
sources and challenge previous reports of a relation between spectrally soft X-ray states 
and the radio emission. 
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1 INTRODUCTION 

The black hole candidate GRS 1915-1-105 was the first Galac¬ 
tic source to display apparent superluminal motions of radio- 
emitting ejecta (Mirabel & Rodriguez 1994; Fender et al. 
1999b; Rodriguez & Mirabel 1999). It is probably the best 
example to date of the strong coupling between the accre¬ 
tion disc and the jet in a black hole system (Pooley & Fender 
1997; Eikenberry et al. 1998, 2000; Mirabel et al. 1998; Bel- 
loni, Migliari & Fender 2000). 

As an X-ray source GRS 1915-1-105 was discovered in 
1992 with WATCH on board GRANAT (Castro-Tirado, 
Brandt & Lund 1992). Early observations, made with this 
instrument and with BATSE/CGRO, showed considerable 
X-ray variability. Since the launch of the Rossi X-ray Timing 
Explorer (RXTE) in 1995 the source has been continuously 
monitored with the All-Sky Monitor (ASM), and frequently 
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observed with the Proportional Gounter Array (PGA) and 
the High Energy X-ray Timing Experiment (HEXTE). From 
these observations the extraordinary X-ray variability dis¬ 
played by GRS 1915-1-105 became clear (see for instance 
Greiner, Morgan & Remillard 1996; Belloni et al. 1997a,b; 
Muno, Morgan & Remillard 1999; Belloni et al. 2000a). Bel¬ 
loni et al. (2000a) found that the entire variability of the 
source is made up of a limited number of distinct variability 
classes; they identified twelve classes that each were observed 
to recur almost identically after intervals of months to even 
years. This behaviour is unique to GRS 1915-1-105, and is 
not observed in any other X-ray binary known to date. 

As in most black hole sources the X-ray spectrum of 
GRS 1915-1-105 can be modelled as a combination of two 
components: a soft (kT^l - 2 keV) disc black body and 
a hard power law extending to >100 keV, generally inter¬ 
preted as representing the radiation from the accretion disc 
and a Comptonizing region (‘corona’), respectively. Belloni 
et al. (2000a) found that the complex X-ray variability of 
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GRS 1915+105 can be reduced to the transitions between 
three basic spectral states which they called A, B and C. The 
spectrally soft states A and B correspond to an observable 
inner accretion disc with different temperatures: in State B 
the inner temperature (~ 2.2 keV) is higher compared to 
State A (~ 1.8 keV). For the spectrally hard State C the 
inner part of the accretion disc is either missing or just un¬ 
observable, and the 2-25 keV spectrum is dominated by a 
power law of spectral index F ~1.3 - 2.4. The description of 
the behaviour of GRS 1915+105 as transitions between these 
three spectral states differs from one based on the canonical 
black hole states (for a review on the canonical black hole 
states, see for instance Tanaka & Lewin 1995 and van der 
Klis 1995), as it is based on the positions in the X-ray colour- 
colour diagram, and not on the combined spectral and tim¬ 
ing characteristics. There are, however, indications that the 
States A, B and C are related to the High State (HS), Very 
High state (VHS) and the Low State (LS) respectively (Bel- 
loni et al. 2000a and references therein). The description 
of the complex X-ray behaviour of GRS 1915+105 as the 
transitions between three basic spectral states, and the clas¬ 
sification of the light curves into twelve different variability 
classes will be discussed in some more detail below. 

GRS 1915+105 is also amongst the most spectacular ra¬ 
dio sources in our Galaxy, displaying apparent superluminal 
motions (Mirabel & Rodriguez 1994; Fender et al. 1999b; 
Rodriguez & Mirabel 1999) and radio oscillations which are 
clearly related to quasi-periodic dips in the X-ray light curve 
on time-scales of tens of minutes (Pooley & Fender 1997; 
Mirabel et al. 1998). The synchrotron emission of these os¬ 
cillations extends well beyond the radio to the millimeter 
and near-infrared regimes (Fender et al. 1997; Eikenberry et 
al. 1998, 2000; Mirabel et al. 1998; Fender & Pooley 1998, 
2000). The power required to generate these repeated syn¬ 
chrotron events is likely to be a very significant fraction of 
the entire accretion energy of the system (Fender & Pooley 
2000, Meier 2001). There is evidence that the power into 
the jet (assuming that a jet is the origin of the radio os¬ 
cillations) may be anticorrelated with the accretion rate as 
inferred from X-ray spectral fits (Belloni et al. 2000b). 

Belloni et al. (1997a,b) suggested that the X-ray dips 
correspond to the disappearance of the inner (< 100 km) 
of the accretion disc. The relation of the radio oscillations 
to these dips led to the assertion that at least some of the 
inner disc which had ‘disappeared’ had in fact been acceler¬ 
ated and ejected from the system, and had not simply fallen 
into the black hole. Based on the similarity between oscil¬ 
lation decays at different wavelengths. Fender et al. (1997) 
came to the conclusion that material is being ejected in the 
form of jets where the synchrotron emitting ejecta are pre¬ 
dominantly subjected to adiabatic expansion losses. What¬ 
ever the interpretation, the coupling between radio-mm- 
infrared emission which is believed to arise in the jet, and 
the X-ray emission from the disc and corona, is very clear, 
and GRS 1915+105 is our most accessible route to date to 
an understanding of the connection between accretion and 
ejection around a black hole. 

In this paper we cross-correlate the entire database of 
radio observations with the Ryle Telescope (RT) at 15 GHz 
over the period 1996 - 1999, with pointed X-ray observa¬ 
tions with the RXTE. In addition we use some data from 
the Green Bank Interferometer (GBI) variable source mon¬ 


itoring program. The next two sections deal with the obser¬ 
vation analysis and the X-ray : radio overlap. In Section ^ we 
present our results by discussing a selection of representa¬ 
tive observations in detail. In Sections and ^the empirical 
characteristics of the X-ray : radio correlation are discussed, 
and Section deals with the profiles of the radio oscillation 
events. In Section we discus a possible scenario which ex¬ 
plains the correlated X-ray : radio behaviour, and finally in 
Section we present our conclusions. 

2 OBSERVATIONS 

2.1 RXTE observations 

2.1.1 Light curve and color-color diagram extraction 

We present 101 X-ray Proportional Gounter Array (PCA) 
observations which overlap with the Ryle Telescope (RT) 
radio observations, made between 1996 June 12 and 1999 
December 13 (see Table |^. Each observation consists of one 
or more intervals (observation intervals) which are separated 
by earth occultations. Following Belloni et al. (2000a), for 
each observation interval we produced three light curves (1 
second time resolution), one in each of the following PCA 
channel bands: a: 0 - 13 (2 - 5 keV), 6: 14 - 35 (5 - 13 keV) 
and c: 36 - 255 (13 - 60 keV). From these light curves we 
calculate two X-ray colours: 

HRi = ( 1 ) 

HR2 = ( 2 ) 

and the total count rate (=a-l-b-i-c). The two colours are de¬ 
fined in this way to ensure linearity: any linear combination 
of two spectral models lies on a straight line connecting their 
locations in the colour-colour diagram (CD; HRi vs. HR 2 ). 
For the background we subtracted a constant level of 10, 
20, and 100 counts in the a, b and c bands respectively, as 
determined from the analysis of typical background spectra. 
Notice that because our data set consists of observations 
from different RXTE gain epochs, we have to consider the 
gain changes. Changes in the PCA gain will have an effect 
on the two X-ray colours resulting in an overall shift in the 
CD. Therefore, gain changes prevent a direct comparison 
between colours from different epochs. However, as the gain 
changes occur on timescales much longer than the typical 
timescales of one observation, we assume here the influence 
on the X-ray colours is negligible. 

2.1.2 Describing the X-ray variability 

Following Belloni et al. (2000a) we analyse the 101 X-ray ob¬ 
servations overlapping with our radio data set, and find that 
most can be assigned to one of their twelve different variabil¬ 
ity classes (see Table ^) . This is shown in Table ^ where we 
give for each observation the MJD, the RXTE observation 
ID-number, the X-ray : radio overlap time in minutes, the X- 
ray class, the length of the short and loM State C intervals, 
and a summary of the radio behaviouijj. State C intervals 

* Observations from RXTE gain epochs 1 and 2 belonging to 
class X were subdivided into four subclasses yi, X 2 , X3 and X4 
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Figure 1. A representative observation in epoch 4 belonging to 
class LJ (observation ID is P-27-00 (#2)). In the top panel we 
show the CD and in lower pan el the light curve. The colours are 
defined as explained in section S- The dips, during which the 
source moves to the lower right in the CD, correspond to state C. 
Outside the dips the source makes rapid transitions between the 
states A and B, reminiscent of the A - B transitions observed in 
class 7 observations (see figure 2 Belloni et al. 2000a). As in all 
the variability classes defined by Belloni et al. (2000a) the source 
makes the transitions in the order A - B - C - A; there are no 
direct transitions from state C to B, but the source always first 
goes through state A (Belloni et al. 2000a). 

are referred to as long or short depending on whether they 
are longer or respectively shorter than 100s. The distinction 
between different classes of variability is made based on the 
light curve and the corresponding CDs; in each variability 
class the pattern in the X-ray light curves and in the CDs 
are the same, while observations from different classes show 
a different pattern (see also figure 2 in Belloni et al. 2000a). 

Most of the observations from our data sample could be 
classified as described above. There were, however, twelve 
observations (see Table and Fig. |^) which showed X-ray 
behaviour that could not be classified into one of the known 
twelve variability classes. Because all the light curves and 
the CDs showed a similar pattern, that did not show any 
resemblance to previously classified observations (compare 
Fig. 0 with figure 2 in Belloni et al. 2000a) we assigned 
these observations to an additional variability class: class ui. 
As already stated by Belloni et al. (2000a), the classifica¬ 
tion of all the observations into variability classes is at this 
stage mostly a tool, allowing one to group similar observa¬ 
tions together and thereby making the complex behaviour 
of GRS 1915-1-105 manageable. 

The description of the observations, in terms of the ba¬ 
sic spectral States A, B and C, is now a relatively straightfor¬ 
ward exercise as the decomposition of observation intervals 

based on the position in the ASM light curve (see Belloni et 
al. 2000a). This subdivision was not made for observations from 
epochs 3 and 4. 


in each of the twelve variability classes into the three states 
was already given by Belloni et al. (2000a). Furthermore 
we can make use of the characteristics of the three spectral 
states as given by Belloni et al. (2000a): 

• State C: low count rate, hard location in the CD 
(HRi-1, HR2>0.1); 

• State A: sharp dip in the light curve, soft location in 
the CD (HRi<l.l, HR2<0.1 ); 

• State B: high count rate, soft location in the CD 
(HRi>l, HR2<0.1). 

Of course, due to the gain changes, the numbers given are 
only valid for gain epochs 1 and 2. The observations from 
epoch 3 used in this analysis consisted only of class x a^nd p 
(see Table and the states could very easily be identified 
from the light curves. In the case of epoch 4 the three states 
still correspond to distinct positions in the CD (although 
these positions are different compared to epochs 1, 2 and 
3), and we found: State A: HRi<0.5, HR2<0.06; State B: 
HRi>0.5, HR2<0.06; State C: HRi~0.4, HR2>0.04. 

2.1.3 X-ray energy spectra 

We extracted energy spectra with an integration time of 16 s 
from Standard2 data in the energy range 3-25 keV (all 
Proportional Counter Units where used when possible, and 
all layers where added). From each spectrum, we subtracted 
the background estimated with Pcabackest 2.1e. Detector 
response matrices were produced using Pcarsp 2.43, and no 
correction for the deadtime was made, as we expect this ef¬ 
fect to be negligible for the purpose of this analysis. We per¬ 
formed fits using the ‘standard’ model for black hole candi¬ 
dates, consisting of the superposition of a multi-temperature 
disc blackbody and a powerlaw component (see Tanaka & 
Lewin 1995). Correction for interstellar absorption and an 
additional Gaussian emission line to take into account an ex¬ 
cess around 6.4 keV were also included. A systematic error 
of 1 per cent was added. 

Belloni et al (2000a) found evidence for a variable Nh 
during the class x intervals (State C only, see Belloni et al. 
2000a) in GRS 1915-1-105. Since leaving the interstellar col¬ 
umn density Nh free to vary during the fits would affect the 
values of the inner radius and make a comparison between 
observations more complicated, we decided to fix it. From 
an analysis of class x observations we found an best fit value 
of 7x10^^ cm"^ 

For the observations presented Section ^ we isolated 
and analysed only the hard, low-flux intervals corresponding 
to State C (see Belloni et al. 2000a). This is because, with 
16 s resolution, we are not able to track the fast variations 
in the remaining parts of the observations. For observation 
K-42-00 there is no State C, but we were able to analyse 
the whole observation. This is also the case for observation 
M-02-00, which consists of State A and C. The observations 
P-30-00/01/02/03 does not consist of any State C, but have 
rapid State A dips. Here again, the fast variations prevent 
reasonable spectral fitting. 

Finally, observation P-16-00/01/02 consists of State 
C only. In this case the energy spectrum is dominated 
by a hard power law component, and the soft (multi¬ 
temperature) disc black body component is not well con¬ 
strained during the 16s spectral fits. This prevented reason- 
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Figure 2. The ~ 4-year radio and X-ray flux history of GRS 1915+105, using daily-averaged data. The top panel indicates the periods 
of overlapping RT and RXTE PCA observations which are presented in this paper. The other panels show the total RXTE ASM count 
rate in the 2-12 keV band, the ASM HR2 colour (5 — 12/3 - 5 keV), radio flux density at 15 GHz as measured with the RT, and the 
2.3 - 8.3 GHz radio spectral index as measured by the GBI. For the spectral index we only show detections at more than 5fT at both 
GBI wavelengths. Also indicated are the epochs of the MERLIN observations of Fender et al. (1999b) and the VLBA observations of 
Dhawan et al. (2000). 


able spectral fitting. However, in order to get an estimate 
of the inner radius we performed fits of each of the three 
individual observation intervals by adding the data of each 
of those intervals together. 

From the spectral hts, performed as explained above, 
the inner radius of the accretion disc is derived from the disc 
blackbody normalization assuming a distance of 11 kpc and 
a disc inclination angle of 66° to the line of sight (Fender et 
al. 1999b). The results are presented in Section Note that 
values of the inner radius are most likely underestimated 
(see Merloni, Fabian, Ross 2000) and may also depend on 
the choice of Nh (as explained above). Therefore, the values 
of the inner radius are only used here as an indication. 

2.2 Ryle Telescope observations 

GRS 1915+105 is monitored regularly at 15 GHz with the 
Ryle Telescope. In bright (more than a few mjy) states data 
can be obtained with up to 8-second time resolution. In this 
paper we use one-day averaged (Fig. ^ and 32-sec radio data 
(Fig. The rms noise on a single 32-sec sample is typically 
5-6 mJy. For further details see Pooley & Fender (1997). 

2.3 Green Bank Interferometer observations 

The Green Bank Interferometer (GBI) monitored variable 
radio sources approximately daily simultaneously at 2.3 and 

8.3 GHz. During periods when GRS 1915+105 is bright 


(^ 10 mJy) it is able to provide high (> 5cr) simal-to-noise 
spectral information, which is utilized in Fig. g. For more 
details see e.g. Waltman et al. (1994). 


3 CORRELATING OBSERVATIONS : THE 

OVERALL PICTURE 

We have cross-correlated the RT and RXTE PGA observ¬ 
ing logs to find simultaneous data, for all observations up 
to 1999 Dec 31. We found 101 overlapping observations, to¬ 
talling 9170 minutes (550 ksec; Table IlJ). It is clear that the 
density of overlapping observations became greater more re¬ 
cently as coordination between RT and RXTE improved. 

Before presenting individual high time-resolution obser¬ 
vations, it is worth looking at the remarkable overall pattern 
of radio and X-ray emission from GRS 1915+105 over the 
past ~ 4 years. This is presented in Fig. 1^, where soft X- 
ray monitoring with the RXTE ASM, the ASM HR2 X-ray 
hardness ratio, RT radio monitoring at 15 GHz, and 2.3 - 
8.3 GHz radio spectral index (from GBI) are plotted. There 
are clearly 6 or 7 major radio flaring events with flares up 
to a few hundred mJy. Most of these radio flares seem to 
be related to a low ASM count rate interval, in combination 
with a high and reasonable stable ASM HR2 level. However, 
not every low ASM count rate, high ASM HR2 interval can 
be related to a radio flare. During the radio flares the radio 
spectral index {a, where Sv oc v°‘) becomes negative, typi- 
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Table 1. List of the 101 simultaneous RXTE PCA - RT observations of GRS 1915+105 up to Dec 31, 1999. The letters I,J,K,L,M,N,0,P,Q 
and R stand for 10408-01, 10258-01, 20402-01, 20187-02, 20186-03, 30402-01, 30703-01, 40703-01, 40403-01, and 40115-01 respectively. 
We give: the observing date (in MJD), the RXTE PCA observation ID, the X-ray variability class as defined by Belloni et al. (2000a) a 
summary of the radio behaviour, and the length of the State C interval (in sec). In column 2 and 4 the “/” is used to make a distinction 
between a sequence of observation IDs. In column 4 a is used to make a distinction between several observation intervals (also referred 
to as “orbits”) for the same observation ID. The radio fluxes, given in column 5, are measured at 15 GHz and when present we give the 
Period (P) of the radio oscillation events. In the last column we define a long State C as an interval during which, on a time-scale of 
~100 seconds, no fast transitions are observed to other states (A and B). The remark “entirely C” refers to observations which consist 
of State C only and which are always longer than 100 seconds. The symbol was used when there are no long and/or short State Cs. 
The observations presented in Fig. W are shown in bold. Note that due to the delay between the radio and the X-ray emission, and the 
relatively short overlaps between the two bands, it is possible that not all the long State C intervals can be associated to radio oscillation 
events just by looking at the table. 


Date 

(MJD) 

RXTE PCA 
OBS-ID 

Overlap 

(min) 

X-ray 

class 

Radio behaviour 
at 15 GHz (RT) 

Length C intervals (s) 

Long Short 

50246 

1-14-01/02 

37 

s 

Undetected (< 1.2 mJy) 

- 

- 

50259 

1-17-03 

109 

<5 

Undetected (< 0.9 mJy) 

- 

- 

50333 

1-33-00 

305 

X4 

Rapid decay 120 —> 20 mJy 

entirely C 

- 

50351 

J-10-00 

148 

h 

Weak (4.6 ± 0.4 mjy) 

- 

34 

50371 

1-41-00 

314 

u 

Weak (2.8 ± 0.3 mjy) 

715 

~100 

50379 

1-43-00 

156 

Xi 

Weak (5.8 ± 0.4 mjy) 

entirely C 

- 

50381 

1-44-00 

206 

V 

Oscillations (P ~ 17 min, AS ~ 60 mjy) 

1100-1900 

~100 

50385 

1-45-00 

50 

X4 

Weak (6.7 ± 0.4 mjy) 

entirely C 

- 

50401 

K-02-02 

154 

X4, V “ 

Weak rise 5 —30 mJy 

220 

- 

50448 

K-09-00 

59 

X2 

Weak (5.6 ± 0.4 mjy) 

entirely C 

- 

50576 

L-01-00 

37 

CL 

Slow modulation 10—30 mJy 

350-2500 


50596 

K-30-01 

50 

X4, OL 

Slow modulation 10—50 mJy 

1950 (x),600 (a) 

~50 

50597 

K-30-02 

2 

CL 

Noisy Sz. variable : 5-35 mJy 

492 

^^80 

50653 

K-39-00 

40 

7 

Weak (1.7 ±0.5) 

- 

- 

50681 

K-42-00 

26 

<5 

Undetected (< 1.2 mJy) 

- 

- 

50698 

K-45-00/01 

440 

0 

Erratic oscillations (P 40, AS" ~ 70 mJy) 

90-540 

~20 - 50 

50706 

M-02-00/01 

540 

e 

Weakly variable 5—30 mJy 

450-750 

- 

50707 

M-02-03/04 

526 

e 

Weak decay 15 ^ 5 mJy 

630-4650 

- 

50708 

M-02-05 

563 

X 

Weak (3.2 ± 0.2 mJy) 

entirely C 

- 

50709 

M-02-06 

101 

X 

Single 50 mJy flare 

entirely C 

- 

50720 

K-48-00 

127 

X4 

Weak (3.8 ± 0.6 mjy) 

entirely C 

- 

50724 

L-03-00 

39 

X 

Bright ~ 100 mjy 

entirely C 

- 

50746 

K-52-00 

240 

X3 

Steady at 60-70 mJy 

entirely C 

- 

50751 

K-52-01/02 

68 

0 

Oscillations (P 22, AS" ~ 90 mJy) 

535-758 

~40-100 

50923 

N-11-00 

141 

X 

Variable ~ 50 mJy 

entirely C 

- 

50945 

N-12-01/02/03 

32 

X 

Variable ~ 50 mJy 

entirely C 

- 

51027 

0-27-00 

49 

p 

Weak (3.0 ± 0.3 mjy) 

- 

'^80 

51036 

0-28-00/01 

159 

p 

Weak (3.3 ± 0.3 mjy) 

- 

'^80 

51040 

0-29-00 

43 

p 

Weak (3.4 ± 0.2 mjy) 

- 

'^80 

51056 

0-31-00 

26 

X 

Weak (10.4 ± 0.3 mjy) 

entirely C 

- 

51071 

0-33-00 

75 

X 

Weak (5.0 ± 0.6 mjy) 

entirely C 

- 

51081 

0-35-00 

128 

X 

Weak (7.4 ± 0.4 mjy) 

entirely C 

- 

51089 

-table continues— 

0-36-00 

65 

X 

Weak (7.6 ± 0.4 mjy) 

entirely C 



Orbit numbers 1,2 and 3 belong to class X4, and orbit number 4 belongs to class v. 
Orbit number 1 belongs to class x and orbit number 2 to class a. 


cally < —0.5, so that flux densities at lower frequencies reach 
even higher values (up to ~ 1 Jy at 843 MHz). Higher time 
resolution radio observations are known to reveal a wealth 
of radio behaviour not apparent in Fig. (e.g. Pooley & 
Fender 1997). 


4 SIMULTANEOUS RXTE PCA AND RT 
OBSERVATIONS 

As already mentioned, Table summarizes all the overlap¬ 
ping RXTE PCA and RT observations of GRS 1915+105 


up to 1999 Dec 31. In Table g we group these into the 13 
X-ray ‘classes’ (consisting of the 12 pr eviou sly known classes 
and the new class defined in Section 2.1.2), and their three 
constituent ‘states’ (Belloni et al. 2000a). Careful examina¬ 
tion of Table ^ reveals that we encounter four different types 
of radio behaviour related to the X-ray classes. First of all 
there are times when the source shows a low level radio flux 
(designated as “weak” with a typical level of a few mjy; 
classes 7 , <5, k, jj,, p and oj) with no radio oscillations. Sec¬ 
ondly, there are observations which have a high radio flux 
level (typically between the 30 - 80 mJy) and show radio os- 
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Table 1 — continued 


Date 

(MJD) 

RXTE PCA 
OBS-ID 

Overlap 

(min) 

X-ray 

class 

Radio behaviour 
at 15 GHz (RT) 

Length C intervals (s) 

Long Short 

51235 

P-07-00 

13 

p 

Weak (4.9 ± 0.5 mjy) 

- 

'^80 

51281 

Q-05-02 

16 

p 

Weak (3.0 ± 0.2 mJy) 

- 

'^80 

51288 

P-12-00 

157 

K 

Weak (3.0 ± 0.5 mjy) 

- 

~20-40 

51299 

P-13-00/01 

191 

7 / 1 ^ 

Weak (2.1 ± 0.3 mjy) 

- 

~100 

51306 

P-14-00/01/02 

132 

K 

Weak (2.6 ± 0.3 mjy) 

- 

~20 

51325 

P-15-00/01/02 

133 

X 

Steady decay 30 —> 15 mjy 

entirely C 

- 

51332 

P-16-00/01/02 

78 

X 

Steady ~ 45 mjy 

entirely C 

- 

51335 

P-17-00/01 

170 

X 

Steady decay 35 —> 10 mjy 

entirely C 

- 

51342 

P-18-00 

141 

p 

Oscillations (P ~ 29, AS ~ 100 mjy) 

130-993 

~40-100 

51352 

P-19-00/01 

113 

p 

Irregular oscillations 0—70 mjy 

1120-1460 

- 

51361 

P-20-00/01/02/03 

211 

X 

Decay 70 ^ 40 mjy 

entirely C 

- 

51370 

P-21-01 

118 

X 

Steady ~10-30 mjy 

entirely C 

- 

51379 

P-22-00/01 

52 

p 

Irregular oscillations 0—70 mjy 

404-1400 

~40-150 

51388 

P-23-00/01/02 

243 

p 

Weak (5.1 ± 0.3 mjy) 

- 

'^80 

51394 

P-24-00 

241 

K 

Weak (3.3 ± 0.4 mjy) 

- 

~10-90 

51399 

P-25-00 

228 

K 

Weak (3.3 ± 0.7 mjy) 

- 

~10-90 

51407 

P-26-00 

227 

K 

Weak (2.5 ± 0.6 mjy) 

- 

~10-90 

51413 

P-27-00 

303 


Weak (3.0 ± 0.2 mjy) 

- 

~100 

51426 

P-29-01 

95 

Lj/ 7 ,tiJ 

Weak (1.7 ± 0.4 mjy) 

- 

~100 

51432 

P-30-00/01/02/03 

118 

k/7 = 

Weak (1.8 ± 0.3 mjy) 

- 

~20-40 

51440 

R-07-00 

53 

7 

Weak (4.1 ± 0.5 mjy) 

- 

- 

51447 

P-31-00 

244 

7 

Weak (1.6 ± 0.3 mjy) 

- 

- 

51465 

P-33-01/02 

125 

(5 

Weak (3.2 ± 0.4 mjy) 

- 

- 

51471 

P-34-00/01/02 

51 

<5 

Undetected (< 1.1 mjy) 

- 

- 

51476 

P-35-00/01 

197 

P 

Oscillations (P ~ 22, AS ~ 60 mjy) 

200-1010 

- 

51484 

P-36-02 

84 

9 

Variable 15-40 mjy 

192-640 

- 

51497 

P-38-00/01/02/03 

151 

eix^ 

Steady decay 140 ^ 40 mjy 

411-1920 

- 

51518 

P-41-00/01/02/03 

117 

X 

Variable 35-85 mjy 

entirely C 

- 

51525 

P-42-01/03 

56 

X 

Variable 55-105 mjy 

entirely C 

- 

M P-27-00 

is classified as class cj. 

, and for P-27-01 orbit number 1 is classified as class 7 and orbit number 2 

as class uj. 

W P-29-01 

is classified as class lo 

, and for P-29-02 orbit 1 is classified as class 7 and orbit 2 as 

class io. 


All the orbits belong to class k, except orbit numbers 2 and 3 from P-30-03 these belong to class 7 . 

(f) P-38-00/03 belong to class 9 and P-38-01/02 belong to class x- 


Table 2. For each variability class, we give 

the average 

maximum radio flux an indication of the radio behaviour. The terms “weak”, 

“variable” 

and “steady” refer to the radio flux level, and the symbol “—” is used when there was 

no overlap between the radio and the 

X-rays. For each variability class 
our data set. 

we also give and their constituent states (Belloni et al. 2000a) and the number of times it occured in 


X-ray 

class 

constituent 

states 

Occurence 

Average radio flux 
at 15 GHz (mJy) 

Radio behaviour 

At 15 GHz 

(y 

ABC 

3 

~ 40 

variable but no oscillations 

p 

ABC 

6 

~ 80 

strong oscillations 

7 

A B 

7 

-- 2.0 

weak, no oscillations 

5 

A B 

5 

^ 1.5 

very weak, no oscillations 

0 

A C 

4 

~ 60 

weak or variable, no oscillations 

K 

ABC 

5 

~ 3.0 

weak, no oscillations 

A 

ABC 

0 

- 

- 

P 

ABC 

1 

~ 5.0 

weak no oscillations 

V 

ABC 

2 

~ 30 

weak or strong oscillations 

P 

ABC 

6 

^ 4.0 

weak, no oscillations 

4> 

A 

0 

- 

- 

X 

C 

25 

~ 5-40 

steady, no oscillations 

(jj 

ABC 

5 

~ 2.0 

weak, no oscillations 
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dilations (classes j3 and u). Thirdly, we have a large number 
of observations (classified as class x) of which about 50% 
have a high level of radio flux (~ 40 mjy) while in the other 
50% it is much lower (;$ lOmJy), but show no radio oscil¬ 
lations. And finally there are observations of class a and 9 
which show a relatively high level of radio emission (~ 40 - 
60mJy) with significant variations, but no radio oscillation 
events. At this point in our argument it is still unclear which 
spectral state (A, B or C) is related to the radio emission, 
however, the only times (4 out of 101 observations) when 
the radio emission was undetectable were in class 5 which 
consist of States A and B only. This suggests that State C 
might be relevant for the production of the (high flux) radio 
oscillation events. 

We now describe in detail several of the high time- 
resolution overlapping observations, in order to clarify this 
behaviour. The observations discussed below are a repre¬ 
sentative subset chosen to highlight the key points of the 
X-ray : radio relation. 


X-ray light curve still displays fairly large amplitude vari¬ 
ability (up to 30 per cent in total intensity), but it is only 
between States A and B, with no State C intervals. In the 
model of Belloni et al. (2000a) this simply corresponds to 
changes in the inner temperature of an accretion disc which 
extends nearly down to the innermost stable orbit. This is 
well illustrated by the fits to the inner disc radius, which 
stay very close to ~ 40 km throughout the observation. Note 
that the temperature changes are in turn presumed due to 
changes in the mass accretion rate through the disc. 

Crucially, there is no radio emission detected during 
the whole observation: the radio flux density is very low, 
and below detectable levels (<1.2 mJy, 3cr). This suggests 
that when there is a bright, spectrally soft accretion disc ex¬ 
tending up to the innermost stable orbit and a hard coronal 
component is weak or absent, the radio emission is low and 
no radio oscillation events are produced (despite apparently 
significant variability in the mass accretion rate through the 
disc). 


4.1 MJD50381 (1996 October 25) 

Presented in Fig. Ha, this observation was first published in 
Pooley & Fender (1997), and was the first direct evidence 
for a disc - jet interaction on short (minutes) time-scales 
from an X-ray binary. The format of Fig. ^ is repeated in 
Fig. c, d, e, f, g and is described in detail in the caption; 
from top to bottem we plot: 

• The radio flux density at 15 GHz, as measured by the 
RT, in 32-second time bins. 

• The total RXTE PCA count rate (with a 1 sec time 
resolution), additionally gray-scale coded to distinguish the 
three basic states (A, B & C), with the class indicated. Em¬ 
phasis is placed on Class C (when present). 

• Best fit values of the inner accretion disc radius (with 
a 16 sec time resolution), as described above. 

When comparing Figs. ^ - g, note that the abscissae and 
the ordinate ranges are not always the same. 

For observation MJD50381 (Fig. |^) we observe a radio 
light curve which is completely dominated by oscillations 
with a quasi-period of ^ 40 min. We will refer to these os¬ 
cillations as the radio oscillation events or simply as radio 
oscillations. The simultaneous X-ray observations reveal a 
light curve (classified as class v) dominated by State C dips 
which appear to recur with the same quasi-period as the ra¬ 
dio oscillations (although the sparse X-ray coverage allows 
for the recurrence period to be twice that of the radio events, 
this is ruled out by observations MJD50698 and MJD51342). 
The (hard) State C dips imply, in the simple X-ray model, a 
‘disappearance’ of the inner ^ 100 km of the accretion disc, 
which then slowly refills. This is reflected by the changes 
in the inner radius, presented in the lower panel of Fig. m. 
Thus we have convincing evidence for the repeated appar¬ 
ent disappearance and refill of the inner accretion disc, with 
associated radio events. 


4.2 MJD50681 (1997 August 21) 

Presented in Fig. these observations reveal a very dif¬ 
ferent relation between disc and jet in GRS 1915-1-105. The 


4.3 MJD50698 (1997 September 7) 

The data from these long (440 min total) overlapping ob¬ 
servations are presented in Fig. The overall classifica¬ 
tion of the X-ray behaviour is class /3: the source clearly 
varies between deep State C dips (reminiscent of MJD50381 
(Fig. ^a)) and rapid A - B - C variations. The deep State 
C dips reach a similar maximum inner disc radius of ~ 100 
km. The pattern of radio emission is such that during the 
periods with deep State C dips, there are repeated radio os¬ 
cillations with a quasi-period of around 20 min, which are 
delayed wit h respect to the X-ray emission by ~30 min (see 
Section 6.3). During the rapid A - B - C variability, however. 


the radio emission drops to a low level, typically ^1.9 ± 0.4 
mJy (which is about the same level as during the State A - 
B transitions seen in observation MJD51432, Fig. ^), with 
no clear bright events. Note that, during the radio oscilla¬ 
tions, the count rate of the State C dips (count rate drops 
to ~ 5 kct sec“^) and the amplitude of the radio oscillations 
(~ 40 mJy) are both comparable to MJD50381 (Fig. ^a). 

The preliminary conclusion that can be drawn from this 
simultaneous observation is that it is not just State C which 
is required to generate high flux radio oscillations, but that 
relatively long State C (longer than 100s, see Section ^ in¬ 
tervals are required in order to generate observable radio os¬ 
cillation events. As the long State C intervals are also both 
deeper and further apart than the short State C intervals, 
any of these three properties (or a combination of them) 
might be required for radio oscillation events. 


4.4 MJD51342 (1999 June 13) 

MJD51342 (Fig. |d) displays another sequence of relatively 
long State C dips followed each time by relatively long peri¬ 
ods of A - B - C variability. In this observation both the ra¬ 
dio and the X-ray repetition periods are clearly visible from 
the light curves. For the X-ray dip recurrence time we find 
^ 30 min. The radio emission shows a very regular sequence 
of radio oscillations with the same recurrence quasi-period, 
unambiguously associating one radio event with each dip. 
It is, however, still not clear to which dip each of the radio 
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0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 

MJD-50381 .0 


Figure 3. a - Simultaneous radio (RT) and X-ray (RXTE PCA) observations of GRS 1915-^105 on MJD50381. The top panel shows the 
RT 15 GHz flux densities plotted with 32-sec time resolution. The middle panel shows the total X-ray count rate in the 2-60 keV band 
(1 sec time resolution), broken up into States A, B and C of Belloni et al. (2000a). The darkest points are those of State C, which the 
analysis throughout this paper indicates is the most significant for jet formation. The class (V’) of Belloni et al. (2000a) is also indicated. 
The lower panel shows the variability of the inner accretion disc radius, Kin, derived from spectral fits to the X-ray data. The X-ray 
dips, apparently caused by the ‘disappearance’ of the inner 100 km of the accretion disc, appear to be closely related to the repeated 
radio oscillations. Figs. are presented in a similar way, but notice that gray scales and both the abscissae and the ordinates ranges 

are not always the same. 



0.79 0.795 0.8 0.805 0.81 0.815 

MJD-50681 .0 


Figure 3. b — Data presented as in Fig. ga, with emphasis placed on state A. Note that in this case there are no State C periods, 
and the inner accretion disc appears to extend right down to ~ 40 km throughout the observation; in this case changes in the flux are 
interpreted as being due to changes in the temperature, which in turn reflects the varying mass accretion rate through the disc. In such 
disc-dominated states, no radio oscillation events are observed and the radio emission is low (typically around a few mjy); in this case 
the radio emission is even below detectable levels. 
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Figure 3. c - Data presented as in Fig. ga. In this case we see that the intermittent p e r iods of radio oscillations appear to occur only 
when there are prolonged, deep State C dips (with a '^30 minutes delay, see Section p.3[ ). During periods of rapid A — B— C flickering 
the radio emission drops to lower (although still detectable) levels and shows no clear individual events. Only during the long dips can 
we make spectral fits to the inner accretion disc radius, which is behaving much as it did on MJD50381 (Fig. ^). 


o 



Figure 3. d - Data presented as in Fig. however, notice the different vertical scale in the top and lower panels. In this case we again 
see a clear relation between long State C dips and radio oscillations. Note that in this case we can clearly see both the radio and X-ray 
repetition periods, which clearly identifies one X-ray dip with one radio oscillation. Note also that in this case the changes in the inner 
disc radius are more dramatic compared to MJD50381 (Fig. and MJD50698 (Fig. ^), as reaches up to ~200 km. 
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LlJ 
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0.75 0.8 0.85 0.9 0.95 1 

MJD-51 432.0 


Figure 3. e - Data presented as in Fig. however, notice that the vertical scale is different in this figure and also that the States A 
and B are highlighted in different colours. Although the X-ray flux shows rapidly repeating soft, State A, dips the 15 GHz radio flux is 
very weak. The X-ray variability is classified as belonging to class k except orbit numbers 2 and 3 from the third interval, these belong 
to class 7 . Because the State A dips are shorter than 16s, no spectral fits could be performed. 



0.95 1 1 .05 1.1 1.15 1.2 

MJD-51 331 .0 


Figure 3. f - Data presented as in Fig. however, notice the different vertical scale in the middle and lower panels. The X-ray light 
curve is classified as class x- the source stays in State C during the whole observation, and only little variation is found in the light curve. 
Not only are there no state transitions to and from State C, also the radio light curve does not show radio oscillation events. Instead the 
radio flux level is high (~ 45 mjy) and almost continuous. Because the X-ray energy spectrum is dominated by a strong hard component, 
the soft disc component is not well constrained which prevented the production of reliable spectral fits. In order to get an estimate of the 
inner radius we performed fits of all the three individual observation intervals added together, and found for the inner radius of 240 
km. 
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0.55 0.6 0.65 0.7 0.75 0.8 

MJD-50706.0 


Figure 3. g - Data presented as in Fig. ga, however, notice the different vertical scale in the top and lower panels. This observation 
belongs to class 6, which shows transitions between State C and A. Although the State C intervals are long enough, there are no clear 
radio oscillations found. Notice that the changes in the inner radius are very small compared to the classes and v. the inner radius of 
the accretion disc is almost the same for state C and A. 


events is related. The X-ray spectral fits indicate that the 
changes in the disc inner radius are more dramatic compared 
to MJD50381 (Fig. |a) and MJD50698 (Fig. ||c), increasing 
up to ~ 200 km. 

4.5 MJD51432 (1999 September 11) 

The X-ray behaviour on MJD51432 is classihed as class w/y 
(see Fig. ^). These classes are characterized by soft, State 
A dips in the light curve which are much deeper than the 
State A dips in observation MJD50681 (about a factor 10 
in count rate lower, see Fig. lb) . Despite the relatively large 
X-ray variability, corresponding to the State A - B transi¬ 
tions, the radio emission is low (~1.8 mjy). These observa¬ 
tions, together with the ones shown in Fig. ^ (MJD50681), 
clearly indicate that (soft) State A and B intervals can not 
be related to radio oscillation events, but instead are char¬ 
acterized by a low-level radio emission. 

4.6 MJD51332 (1999 June 3) 

The previous examples have shown either the case when 
there is a low level radio emission and no radio oscillation 
events, or when there is a high radio flux and radio oscilla¬ 
tion events. In Fig. ^ we show an example of an observation 
which has a relatively high level of radio flux, but no radio 
oscillations. The X-ray light curve is classified as class x ^-nd 
shows very little variations and/or state transitions. In fact, 
the source stays in State C for periods of days to months. 
From these observations we conclude that an uninterrupted 
State C interval can produce a high level of radio emission, 
but no radio oscillation events. However, note that there 


are a considerable number of uninterrupted State C inter¬ 
vals with low radio flux (^5-10 mJy), also without radio 
oscillation events (see Table M. 

As explained in Section p.l.j we can only get a rough 
estimate of the inner radius of the accretion disc from the 
spectral fits due to the fact that the soft component is not 
well constrained. However, from spectral hts of the whole 
observation (all intervals added) we hnd for the inner radius 
of the accretion disc ^ 240 km. This suggest that the inner 
radius is much larger during the class y observations than 
during the long (interrupted) State C intervals. 

4.7 MJD50706 (1997 September 15) 

The observation presented in Fig. ^ is classified as class 
9 and shows transitions between State A and C. When we 
compare this observation with the ones presented in Fig. |^, 
c and d we immediately notice two major differences. First of 
all, although the State C intervals in class 6 are long enough 
(^ 100s, they are comparable in length to the State Cs in 
class P and v) and seem to be related to a relatively high 
level of radio flux (see Table only strong variability (with 
no periodicity) is found in the radio flux and clear oscillation 
events such as found in class P and v are absent. Secondly, 
the radius hts do not show the same pattern (fast increase 
in the inner radius, followed by a slower decrease) during 
the State C intervals we found in class P and v. In fact we 
hnd that the inner radius shows no large variations and the 
value we hnd in State C is comparable to that in State A. 

From Table ^ and ^ it is clear that class 9 is the only 
class which shows long State C intervals (like the ones from 
class P and u which are related to the radio oscillation 
events), but no State B intervals. Another difference be- 
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Table 3. The Linear correlation coefficient and the Spearman Rank correlation coefficient for the correlations between the RT radio flux 
and the RXTE ASM flux and HR2 colour, as shown in Fig. A positive value of the correlation coefficient indicates a correlation, and 
a negative value an anticorrelation. For each of the correlations the confidence level (in percent) is given. 


Radio vs. X-ray 

Coeff. 

XI 

Conf. 

Coeff. 

X 2 

Conf. 

Coeff. 

X 3 

Conf. 

Lin. corr. coeff. 

-0.63 

> 99% 

+0.52 

> 99% 

-0.87 

> 99% 

Rank corr. coeff. 

-.044 

~ 98% 

+0.63 

> 99% 

-0.82 

> 99% 

Radio vs. HR2 


XI 


X 2 


X 3 


Coeff. 

Conf. 

Coeff. 

Conf. 

Coeff. 

Conf. 

Lin. corr. coeff. 

+0.47 

~ 99% 

-0.14 

~ 68% 

+0.88 

> 99% 

Rank corr. coeff. 

+0.38 

~ 95% 

-0.15 

~ 79% 

+0.84 

> 99% 


tween class 6 and classes /? and v is that the long State 
C dips in class 9 are near-contignous, with only very little 
time (;$ 3 min) between them, whereas in classes j3 and v 
this is much larger (~ 20 min). So, another possibility is 
that the additional requirement on State C dips to produce 
radio oscillation events is that they are “well-separated”. 


5 CLASS X OBSERVATIONS 

Besides the long State C intervals which are related to the 
radio oscillation events, there are also times when the source 
occupies State C for periods of days to months. Observations 
during these periods are marked “entirely C” in Table and 
are characterized by steady radio and X-ray emission. An 
example of such an observation is shown in Fig. Belloni 
et al. (2000a) classified all these observations, available to 
them at that time, as class x and subdivided them, as a time 
sequence, in the subclasses xi, X 2 , Xs- Below we will discuss 
the three subclasses in more detail. Also a fourth subclass, 
X4, was introduced which comprises all the observations that 
did not fall in one of the three other classes. As this subclass 
does not consists of a time-sequence of state C observations, 
it will not be discussed as a separate class below. 

5.1 The radio behaviour during the x classes 

The intervals xi and X3 are referred to as “plateau” states 
(or as radio-loud, radio-plateau (Muno et al. 2001), or type II 
states (Trudolyubov 2001)) and a handful of these states 
have been observed since 1996 (Pooley & Fender 1997, 
Fender et al. 1999b, Fender 2001). The pattern of these 
“plateaux” seems fairly stable (Pooley & Fender 1997; Han- 
nikainen, Hunstead, Campbell-Wilson 1998; Fender et al. 
1999b, Fig. They are generally preceded by (an) opti¬ 
cally thin radio event(s) (“pre-plateau flare”), settle into the 
quasi-stable state for a few weeks (“plateau jet”), and then 
are followed by further optically thin radio event(s) (‘post¬ 
plateau flare’; see also Fig. |^. During the quasi-stable state 
the radio emission seems to saturate at ~100 mjy (at 15 
GHz; Fig. ^), while the radio flux during pre-plateau and 
the post-plateau flares generally is a factor 3 to 6 higher 
(see top panel Fig. ^ and Mirabel & Rodriguez 1995, Han- 
nikainen et al. 1998 and Fender et al. 1999b). It is only for 
X3 that the quasi-continuous ejections and the discrete ejec¬ 
tions corresponding to the quasi-continuous radio jet and 


the radio flares respectively are directly observed (Dhawan 
et al. 2000). 

During class X2 the radio flux is typically a factor 10 
lower compared to xi a-nd X3 (see Fig. and no major 
ejections of relativistic plasma have been reported. Hence 
the X2 interval is often referred to as radio-quiet (Muno et 
al. 2001; or as type I (Trudolyubov 2001)). 

Simultaneous infrared coverage in 1996 (Bandyopad- 
hyay et al. 1998) indicates that the flat (q;~ 0, at GHz fre¬ 
quencies) radio spectrum during the plateau states may ex¬ 
tend to near-infrared wavelengths (Fender 2001). In Fender 
et al. (2000) the fiat-spectrum emission of GRS 1915-1-105 
during the plateau states is compared to that of Cyg X-1 
and Cyg X-3, and it is argued that such spectra arise in 
a conical, quasi-continuous, partially self-absorbed outflow 
such as originally envisaged (in an idealized case) by Bland- 
ford & Konigl (1979) and developed by others, including 
application to X-ray binaries (e.g. Reynolds 1982; Hjellm- 
ing & Johnston 1988; Falcke & Biermann 1996, 1999). In 
this case, the plateau states correspond to the formation 
of a luminous quasi-steady outflow from GRS 1915-1-105, a 
scaled-up version of the jet from Cyg X-1 in the Low/Hard 
X-ray state (Stirling et al. 2001). 

For class X2 we combined the GBI and Ryle telescope 
averaged data points to obtain the first reported radio spec¬ 
trum during X2. The flux at 2.3 GHz showed large variations 
as a function of the Hour Angle (HA), while the 8.3 GHz 
flux stayed almost constant. The larger field of view at 2.3 
GHz may result in a significant contribution to confusion 
from nearby sources, and this may result in the fluctuations 
noted. We therefore use the data at 8.3 and 15 GHz. The av¬ 
erage radio spectrum is plotted in Fig. P, and for the spectral 
index (between 8.3 - 15 GHz) we find: a=—0.2 ± 0.1. This 
means that although during X2 no radio jets or radio os¬ 
cillations were observed and the radio flux was significantly 
lower compared to classes xi and X3, the (high frequency) 
radio spectrum is also found to be approximately fiat. 


5.2 X-ray and radio correlations 

In Fig. 1^ we compare daily-averaged radio and X-ray data 
for the prolonged State C periods Xi) X2 & Xs- In Table ^ we 
show both the Linear and the Spearman Rank correlation 
coefficients (Press et al. 1992) and their confidence level for 
all the relations shown in Fig. p. From this we can conclude 
that for the plateau states xi nnd xs: 
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Figure 4. A schematic view of the radio oscillation events, the pre- and post-plateau flares and the plateau jets in GRS 1915-1-105. 
In the top panel we show the radio (RT and GBI) and the X-ray (RXTE ASM) light curves from MJD50700 to MJD50780. We have 
indicated the positions of the flares, the jet and also of the start of the radio oscillation events (which start just after the beginning of the 
post-plateau flare). Also indicated are the times the source was observed with MERLIN (see also Fender et al. 1999b). In the lower panel 
we show schematic views of the ejection events we encountered (the numbers correspond to the positions in the light curve indicated 
in the top panel): 1) the pre- and post-plateau flares; 2) the plateau jet, with quasi-continuous ejections and no single blobs observed; 
3) the radio oscillation events, a sequence (in the direction of the arrow) is shown during which a new blob is ejected. The figures are 
not drawn to scale and estimates for the sizes of the jets and the corona are indicated. Rs is the Schwarzschild radius, and we assumed 
a black hole mass of 10 Mq. For the pre- and post-plateau flares the blobs are observed up to distances of 3000 AU from the core 
(Fender et al. 1999b, Dhawan et al. 2000). In the case of the plateau jet the distance to the optically thin point (t=1) at 15 GHz is 
given (see also Dhawan et al. 2000), and for the radio oscillation events a typical distance (upper limit) between two blobs is estimated 
from the period of the oscillations ('^ 30 min) and assuming the blobs are moving with the speed of light. Notice that during the radio 
oscillation event sequence each new ejection starts before the previous one has completely faded away, as can be found from the light 
curves (Figs. c, and d). 
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Figure 5. From left to right: the 15 GHz radio flux density as a function of RXTE ASM 2-12 keV flux, the 15 GHz radio flux density 
as a function of the hardness ratio (HR2 = RXTE ASM 2-12 keV flux as a function of the hardness ratio, using 

daily-averaged data for the long State C intervals yi, X2 and X3 (Belloni et al. 2000a). Note that for xi X3 there is anticorrelation 
between radio flux density and total X-ray flux, and a correlation with hardness. For X2 the situation is reversed: the radio emission is 
correlated with X-ray flux and anticorrelated with hardness. For all the three x intervals there seems to be an anticorrelation between 
the X-ray flux and the hardness ratio: when the source becomes brighter, it gets softer. 


• the radio flux density is anticorrelated with the soft X- 
ray flux; 

• the radio flux density is correlated with the ASM HR 2 
colour (for definition see caption Fig. ^). 

And for the state X 2 we find: 

• the radio flux density is correlated with the soft X-ray 
flux; 

• the radio flux density may be anticorrelated with the 
ASM HR 2 colour. 

So while the plateau states xi s^nd X 3 show very similar be¬ 
haviour with respect to the radio and X-ray correlation, the 
X 2 interval shows a different behaviour. Note that this divi¬ 
sion is also reflected in the X-ray behaviour: X 2 has a lower 
count rate compared to xi and X 3 (Fig- 1^- Also Belloni et 
al. (2000a) found that X 2 occupies a lower branch in the CD 
compared to xi and xs, a phenomenon they explained as 
related to a difference in N_tf and/or a difference in the high 
energy cutoff. 


O 
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Figure 6. The radio - infrared spectra for the plateau jet, the 
X2 interval and the radio oscillation events. The X2 spectrum is 
the first reported radio spectrum for this interval. The sp ectr um 
is flat, with a spectral index of q:=— 0.2 zb 0.1. See Section ^.l| for 
more details. The spectrum of the plateau interval is taken from 
Fender (2001), and the spectrum of the radio oscillation events is 
taken from Fender Pooley (2000). 


6 EMPIRICAL FEATURES OF THE DISC - 
JET CONNECTION 

6.1 State C intervals 

In order to check that indeed we need a long, single State C 
interval to produce a high level of radio emission, and not 
just ‘a lot of State C’ we calculated the time the source was 
in State C during the intervals when no long State Cs where 
present for observation MJD50698 (Fig. ic). We compared 
this with the time the source shows long State C intervals 
(for the same observation). We find that for the intervals 
around MJD50698.67 and MJD50698.81, which are associ- 


atedQ with a high radio flux (and radio oscillation events), 
the source shows long State Cs for 45 percent and 27 percent 
of the total time, respectively. On the other hand, the inter¬ 
vals around MJD50698.74, MJD50698.94, MJD50698.98 and 
MJD50699.0, which are related to a low-level radio emission, 
show short State C for 41 percent, 30 percent, 29 percent and 
31 percent of the total time, respectively. This clearly shows 
that the presence alone of State C is not enough: there are at 

t Here we used the fact that the radio e rniss ion is delayed with 
respect to the X-ray emission, see Section |6.3| 
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least four observation intervals which lack the long State Cs, 
the high radio emission, and the radio oscillation events, but 
spend more or a comparable fraction of their time in State 
C compared to the observation intervals which do show long 
State Cs. It also shows that the long State C intervals have 
a disproportional higher level of radio flux compared to the 
short State C intervals. 

Is it only the length of the State C intervals that de¬ 
termines if the radio oscillation events are produced? From 
Section and Table ^ it is clear that classes 6 and a do 
show long State C intervals and a relatively high level of ra¬ 
dio flux, but no clear radio oscillation events. Unfortunately, 
due to the poor quality of the X-ray : radio overlap in the 
class a observations, we were not able to say anything about 
the behaviour of the source in this class. For class 0 we do 
have good overlaps, and Fig. ^g shows an example. From 
Section 4.7 and the above, it seems that for the high flux 
radio oscillation events to be observed, two conditions have 
to be met: 


• the State C intervals have to be longer than ~ 100s; 

• the State C intervals have to show a repeating pattern 
in which they are either separated by State B intervals, or 
“well-separated” in time, or both; 


As can be seen from the Figs c and d the long State C 
intervals are also the ones which have the lowest X-ray count 
rate. Therefore, the first statement that long State C inter¬ 
vals are related to the radio oscillation events also implies 
that low-count rate State C intervals are related to the radio 
oscillation events. It would take observations of short State 
C intervals of low count rate to decide between these two 
possibilities, but such events have never been observed yet. 
Regarding the second statement, it is important to note that 
it is not clear at this point why we do see radio oscillation 
events in classes (3 and v and not in class 8: it can be due 
to the presence (or absence) of the soft State B intervals, 
or to the larger (or smaller) separation between the State 
C intervals, or both. Although we feel that the presen ce o f 
State B intervals can be very important (see Section 8.2), 
the most simple solution is that the State C intervals need 
to be well-separated for the individual radio events to be 
observed. 

It is interesting to note the difference between State Cs 
of class 8 and the State Cs of the other variability classes. 
Besides the differences already mentioned in Section 4.7, we 
also find class 8 to be different in several other ways. From 
comparing the lightcurve in Fig. ^ with the ones from class 
(3 and u (Fig. |^, c and d) it is clear that not only do the 
State C intervals of class 8 have the highest count rates com¬ 
pared to the State Cs in other classes, but the count rates 
are also higher than State A intervals. Regarding the spec¬ 
tral parameters, the class 8 State Cs are also different. From 
the spectral hts we found an inner radius of the accretion 
disc which is about 3 times smaller compared to the State Cs 
in classes (3 and u (see Figs. c, d and g), while the inner 
temperature is comparable to what is found during State A 
and B (~ 1.5 keV). Together with the fact that the power 
law component is very steep (F ~ 3), we would conclude 
that these State Cs are spectrally softer. However, from the 
best ht parameters we find that on average the luminos¬ 


ity 0 of the power law component (~ 2.2x10®®erg s ®) is 
larger than the luminosity of the disc blackbody component 
(^ 1.9x10®®erg s“^). For the State C intervals in classes (3 
and V on the other hand, the power law luminosiy is smaller 
(^ 0.7x10®®erg s“^) compared to the disc blackbody lumi¬ 
nosity (~ 2.5x10®®erg s“®). Although it is obvious that the 
State C intervals are clearly different in class 8 compared to 
the other classes and we cannot rule out the possibility that 
these differences are responsible for the “anomalous” radio 
behaviour, more analysis is necessary to say something more 
firm. 


6.2 State A and B intervals 

One can ask the question whether it is only long State C in¬ 
tervals which produce radio oscillation events, or if perhaps 
long state A and B intervals do the same; i.e. is it perhaps 
only the length of the interval which determines the produc¬ 
tion of the radio oscillation events? The strongest evidence 
against long State B intervals being responsible for the pro¬ 
duction of radio oscillation events is shown in Fig. In this 
case we clearly have long (up to ^ 30 min) state B intervals, 
but the radio emission is below detectable levels. The only 
cases where we have long state A intervals are during the 
observations presented in Figs. d and g. These observa¬ 
tions are classified as class jS and 8 and in fact it is only 
for observations in these classes and class (j) (which consists 
of State A only, but for which we do not have any overlap¬ 
ping radio observations) that the source shows long State 
A intervals. In classes (3 and 8 each long State C interval is 
followed by a (slightly shorter, ~ 80 - 100s) State A inter¬ 
val. Hence, one might relate these long State A intervals to 
the radio oscillation events. However, when we compare the 
observations presented in Figs. ^ and d with the class u ob¬ 
servation shown in Fig. ^a, which does not have long State A 
intervals, we see the same pattern of radio oscillation events 
(although the pattern in Fig. ^g is not so clear, we do be¬ 
lieve the State C intervals are responsible for the variations 
in the radio flux; see above). We therefore conclude that 
only (>100s) long, well-separated State C intervals produce 
radio oscillation events. 

What is the difference in the radio behaviour when the 
X-ray light curve consists of short State C dips and when 
it consists of short State A dips? From Fig. and Table ^ 
It is clear that in both cases the radio emission is low. We 
find that the average flux level corresponding to intervals 
showing short State A dips reaches 2.0 ± 0.10 mjy while for 
the intervals showing short State C intervals this is 3.3 ± 
0.10 mJy (numbers calculated from Table |^. From this we 
conclude that, although they are clearly not related to the 
radio oscillation events, short State A intervals and short 
State C intervals show a very similar radio behaviour. In 
fact State A, B and short State C intervals are all related to 
low-level radio emission. For some observations we see a slow 
decay in the low-level radio flux and therefore, at this point 
we cannot rule out the possibility that the radio emission is 
the decay from a large radio event prior to the observation 
in question, as these radio events can be very short (^ 24h). 




Luminosities are calculated assuming a distance of 11 kpc. 
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6.3 The one — to — one relation; estimating the 
time delay 

It is clear from the figures that the X-ray and radio events 
do not occur simultaneously. In Fig. the radio oscilla¬ 
tions switch ‘on’ and ‘off’, allowing us to relate each State 
C interval with one radio oscillation event. This is indicated 
in Fig. ^ where we show a close-up of the X-ray and the 
radio light curves of MJD50698. Despite the sparse cover¬ 
age in the X-ray due to earth occultations, in view of the 
known characteristic repeatable X-ray behaviour we believe 
there is convincing evidence for a one-to-one relation (as in¬ 
dicated in Fig. ^ between the dips and the radio oscillation 
events. Moreover, we can also associate each X-ray dip with 
one of the radio oscillation events occurring after this X-ray 
dip. This allows us to relate each of the State C intervals 
in the X-ray light curves to one radio oscillation event in 
the other observations presented here. A strong indication 
that the one-to-one relation is indeed a intrinsic feature of 
the source and not just a coincidence, is the fact that the 
periods found in the radio oscillation events are the same as 
the periods of the long hard State C dips in the X-ray light 
curve: we find periods around 30 minutes in good agreement 
with the range of 20 - 40 min found by Pooley & Fender 
(1997). Furthermore, the one-to-one relation allows us to 
“fill in the gaps” in the X-ray data which are due to earth 
occultations. If we assume the source shows a continuous 
behaviour in each class and we can safely fill in the missing 
gaps, this does not lead to any inconsistencies. An excep¬ 
tion to this seems to be class 9, however, we do feel that 
also in this case each State C interval is responsible for an 
ejection event and the subsequent radio emission, although 
we can not distingui sh b etween these individual events (as 
explained in Section 6.1). 

The one-to-one relation allows us to estimate the time 
delay between the radio and the X-ray emission. From Fig. ^ 
we find a delay of ~30 min; similar delays (14-30 min) were 
also found in the other observations presented here. These 
delays are in agreement with theory: the self absorption in 
the optically thick region of the jet causes the radio emission 
to be delayed compared to higher frequencies, e.g. IR. Note 
that we measure this time delay from the end of the State C 
to the peak of the first subsequent radio event. These refer¬ 
ence points were chosen because these are easily determined 
from the figures. This way we can only get a lower limit 
of the time delay, as the real starting point of the ejection 
(which does not have to be the end of State C) can not be 
determined. Note that measuring the time delay from the 
beginning of State C to the peak of the subsequent radio 
event would make the delay longer by ~ 10-15 min, which 
is the average length of the State C intervals. 


7 OSCILLATION EVENTS - DELAYS AND 
PROFILES 


By comparing simultaneous radio and IR observations 
Fender & Pooley (1998) found that most likely the IR pre¬ 
cedes the radio by ~30 min, which is about the same time 
dela y as found for the X-ray and the radio emission (Sec¬ 
tion O). This infers that the delay between the X-ray and 
the IR is small, a result already found by Mirabel et al. 
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Figure 7. Close-up of radio and X-ray light curves on MJD50698, 
indicating our best guess at the one-to-one correspondence be¬ 
tween State C dips and radio oscillation events (see also Fig. pp). 


(1998). Mirabel et al. (1998) fit the frequency-dependent 
time delays from the near-infrared to radio wavelengths with 
a simple van der Laan (1966) model. That model predicts 
the time delay to peak emission at frequency v with a func¬ 
tion t^ oc (depending on the electron distribution in¬ 

dex) where tv is the time delay. In an idealized conical out¬ 
flow, the characteristic size scale varies with frequency as 

(e.g. Blandford & Konigl 1979). If the flow is already 
accelerated to its terminal velocity by the time we observe 
the oscillation events then the time delay tv oc . For 
nearly all conceivable models in which optical depth effects 
cause the frequency-dependent delays in the oscillations, the 
infrared emission is likely to arise within seconds of what¬ 
ever is the disc trigger event (due to lower optical depth at 
higher frequencies). Therefore, directly associating jet for¬ 
mation with some particular part of the X-ray dip may be 
more profitably pursued at infrared wavelengths. However, 
the following should be borne in mind with respect to the 
emission profiles observed at different wavelengths: 

• The similarity of the oscillation profiles at radio and 
infrared wavelengths, and the overall flat radio-infrared 
spectrum, both indicate that the self-similar phase of syn¬ 
chrotron source evolution has already been reached by the 
time the infrared emission is observed. 

• This implies that optical depth effects do not play a 
role in determining the observed profiles of the oscillation 
events, and that all IR-mm-radio events are observed from 
material already moving at, or near to, the terminal velocity 
of the outflow. 

If the rise phase of the oscillations is not an optical 
depth effect, it probably either reflects: 

(1) The period of Injection of relativistic particles Into the 
outflow from the dlsc/corona 

(li) The ‘width’ of a shock which is passing down the jet 
accelerating particles in-situ (e.g. Kaiser, Sunyaev, Spruit 
2000 ). 

In case (i) we may well be able to relate the rise timescale to 
some timescale in the X-ray light curve; in case (ii) this is not 
so clear, since we do not understand how the shock would 
be generated. However, in both of the above cases (i and ii), 
the lack of optical depth effects in the rise profile indicates 
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that the timescale for particle injection/acceleration is long 
compared to the time required for the flow to pass through 
its T = 1 surface (“photosphere”). 

Moreover, based on the similarities in the oscillation 
profiles at different wavelengths, the decay phase proba¬ 
bly reflects the intrinsic adiabatic decay of the expanding 
cloud which is also “smeared out” due to the relativistic 
Doppler effects. As the decay phase is not related to the X- 
ray behaviour a comparison of the IR-mm-radio decay phase 
with the X-ray light curve is not relevant (regardless of the 
Doppler effects). 

All of this leads us to conclude that: 

• we may well be able to usefully compare the rise time 
of the radio oscillation events with some timescales observed 
in the X-ray light curve; 

• both the rise phase (particle injection and/or acceler¬ 
ation by passing shock) and decay phase (adiabatic expan¬ 
sion losses) are frequency-independent effects, which nicely 
explains why the profiles are identical at all frequencies. 

We note three caveats to the above line of reasoning - 
firstly, at lower radio frequencies the oscillations do appear 
to become a bit less sharp, perhaps indicating the optical 
depth effects are beginning to play a role (see e.g. Mirabel 
et al. 1998); secondly, in at least one model the hard X-rays 
are also beamed (Markoff, Falcke & Fender 2001), although 
not by a very large factor; thirdly, besides the adiabatic ex¬ 
pansion the decay phase may additionally be affected by the 
decay of the injection of material into the jet (which is not 
affected by the Doppler effects). The above line of reasoning 
may well help to constrain the physical dimensions and op¬ 
tical depth structure of the jet in GSR 1915-1-105 and other 
systems. 

It is worth making a further clarification here. When we 
state ‘ejection’ or ‘injection’ we are talking about a popula¬ 
tion of relativistic particles which we observe via their syn¬ 
chrotron emission. We cannot rule out possibilities such as a 
continuous outflow of particles, which in soft states (A and 
B) are cooled before we can observe the synchrotron emis¬ 
sion. This cooling may occur via the inverse Compton pro¬ 
cess, where the nascent population of relativistic electrons 
would lose their energy before having a chance to propagate 
downstream in the jet to regions from which the opacity 
dropped sufficiently for the synchrotron emission to be ob¬ 
servable. In states A and B the soft X-ray spectrum peaks 
at around ^ 1 keV, and the synchrotron emitting electrons 
probably have Lorentz factors 7 > 100 (e.g. Fender et al. 
1999b) and so the inverse Compton process would scatter 
the photons out of the easily observable X-ray band (the 
increase in energy is typically by a factor of order 7 ^). Is the 
photon field sufficient to cool these electrons? A detailed 
calculation of this is beyond the scope of this paper, but a 
simple estimate of the energy density of the photon field 50 
Schwarzschild radii (a somewhat arbitrary guess at where 
the jet might be launched) above a 10 ^® erg s“^ inner disc 
region is of order 10 ^° erg cm“® (assuming a black hole mass 
of IOMq). In such a photon field the lifetime of an electron of 
initial Lorentz factor ~ 100 to inverse Compton losses would 
be 10 “® s, which is not long enough for the electrons to es¬ 
cape to regions of significantly lower photons densities. Note 
that the spectral steepening related to the inverse Compton 
cooling has never been observed yet; this could be related 


to the fact that the transition between jet and cooled states 
is so rapid that we don’t have a chance to catch the onset 
of cooling. 


8 DISCUSSION 

In this section we discuss the results of our analysis. Our 
main concern is with the radio oscillation events, however, 
the radio emission (pre- and post-plateau flares and the 
quasi-continuous plateau jets) observed during the plateaux 
states will also be briefly discussed. In the following when¬ 
ever we refer to the term ‘jet’, this includes the pre- and 
post-plateau flares, the quasi-continuous plateau jets and 
the radio oscillation events, which we assume are all related 
to an outflow. 


8.1 The disc—jet connection 


The analysis of the complete set of overlapping radio and X- 
ray observations of the black hole candidate GRS 1915-1-105, 
presented in the previous sections, has unambiguously 
shown the intimate relation between the radio oscillation 
events and the spectrally hard State C episodes in the X- 
ray time series. As already shown by Belloni et al. (2000a), 
State C is the state in which GRS 1915-1-105 is most similar 
to the canonical black hole Low/Hard state (although it is 
not completely identical). A relation between the Low/Hard 
state in black hole X-ray binaries and their radio emission 
(quasi-continuous jets) has already been established previ¬ 
ously (Fender 2001): in the Low/Hard X-ray states of at least 
three systems (Cyg X-1, IE 1740.7-2942 and GRS 1758-258) 
quasi-continuous radio jets are directly imaged, whereas a 
general feature of all ‘jet-systems’ in the Low/Hard state 
(see Fender 2001, Meier 2001) seems to be a flat or inverted 
synchrotron spectrum extending from radio up to near IR 
or even sometimes optical wavelengths. GRS 1915-1-105 is no 
exception to this rule, as both the flat spectrum (Fig. ^) and 
the plateau jets themselves are (directly) observed (Dhawan 
et al. 2000 ) during the plateau states. 

During the hard state (G) the X-ray spectrum is domi¬ 
nated by a power law component, which is widely accepted 
to originate from inverse Gompton processes in the inner 
geometrically thick accretion flow (such as a corona and/or 
an ADAF), although there are alternatives (e.g. Markoff et 
al. 2001). The empirical link between the radio emission 
and these states, especially the one-to-one relation (see Sec¬ 
tion 3.3), strongly suggests that besides the outer geomet¬ 


rically thin, and inner geometrically thick inflow there is an 
outflow in the form of a jet. Fender et al. (1999a) have sug¬ 
gested that during the Low/Hard state in black hole systems 
the Gomptonizing corona and the radio jet are strongly cou¬ 
pled. Also, based on a study of the spectral and the timing 
properties of GRS 1915-1-105 during the plateau states and 
the radio quiet states (X 2 ), Muno et al. (2001) come to the 
conclusion that a coupling between the corona and the jet 
seems very likely. 

Assuming the jet and the corona are coupled, an injec¬ 
tion of material into the jet from the corona (or generation 
of a shock - in the future we will not repeat this alternative 
but consider it implicit) takes place during State G intervals. 
If this ejection occurred at only one point, for instance as 
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the source moves in or out of State C, we would expect the 
rise time of the radio emission to be shorter compared to the 
case when the ejection takes place over a finite interval in 
time. In fact from the arguments given in Section the fi¬ 
nite rise time indicates that the particle injection cannot be 
a ‘delta function’ but must also have a finite, observable du¬ 
ration. From the observations shown in the previous sections 
we find the rise time of the radio to be broadly comparable 
to the length of the State C intervals (~ 10 - 15min), in¬ 
dicating that a more or less continuous ejection takes place 
during State C. 

From the considerations given above the following sce¬ 
nario emerges for the production of the radio oscillation 
events (see also sequence shown in Fig. ^ panel 3). When 
the source moves into State C the thin disc truncates, as the 
inner radius moves outward (in the case of Fig. c and d, 
this is occurs within 3-5 min). In State C there is a large 
corona, and coupled to it a jet. The outflow through the 
jet comes about through the more or less continuous ejec¬ 
tion of material from the corona. As a direct consequence of 
this, longer State Cs are also radio brighter as more mate¬ 
rial is being ejected. As the source moves out of State C the 
accretion disc moves back in, a process which occurs on a 
longer timescale (~ 9 - 25 min, see Fig. ^a, c and d), even¬ 
tually causing the injection of relativistic particles to cease. 
The current radio data are not really of high enough quality 
to see if there are measurable changes in the injection rate 
(ie. the data seem to be consistent with a more or less con¬ 
stant gradient in the radio light curve). Higher quality data 
may also be able to detect the effects of the decay phase 
emission from both approaching and receding components 
(assuming the ejections are symmetric), which will have dif¬ 
ferent Doppler factors. The result of the X-ray state change 
is that we observe a discrete volume of synchrotron-emitting 
particles moving downstream in the jet. 

For sequences of repeating long (^ 100s) State C inter¬ 
vals we only observe the radio oscillation events when the 
separation between the State C intervals is large enough for 
the individual events to be observed: if the time it takes for 
the material to move up into the jet (from optically thick to 
thin) is longer than the recurrence time of the State C inter¬ 
vals, we can not distinguish between the individual events. 
The length of the State C intervals is not important in that 
respect: the length only determines the strength of the radio 
emission. 


8.2 Long and short State C intervals 


In light of the scenario for the prod uctio n of the radio oscil¬ 
lation events, suggested in Section 8.1, the question arises 


why series of short State C intervals do not also produce 
(high flux) radi o os cillation events? 

In Section | 6 .l| we have shown that for the (high flux) 
radio oscillation events to be observed, two conditions have 
to be met. The short State C intervals are not only too short 
100 s) to produce a high level of radio emission, but also 
the separation between the State C intervals is not enough 
for the individual events to be observed. 

To explain the disproportional high radio flux for the 
long State C intervals compared to the short State C inter¬ 
vals, we suggest that cooling of the material in the jet by 
soft X-rays from the State A and/or B (which are found in 


between the State C intervals; Section might be impor¬ 
tant. In the case of series of short State C intervals the soft 
photons of the disc can catch up with the material in the 
jet before it reaches the r = 1 surface, as the State Cs are 
both short and close together, and hence the radio emission 
is suppressed. 


8.3 The low-level radio emission 


It has become clear that State A, B and short State C in¬ 
tervals are related to a low-level radio emission (see Fig 
The short State C intervals have already been discussed in 
the previous section. But why do we see radio emission in 
States A and B when there is no hard State C present (like 
in the case of class 7 and 5)7 

One possibility is that the low-level radio emission we 
observe is indeed related to the decay of a large event occur 
before the interval in question, as noted in Section ^ An¬ 
other possibility is to realize that the absence of State C in 
the X-ray light curve does not necessarily mean that there 
is no hard component in the X-ra y spectrum. In fact, as al¬ 
ready mentioned in Section 2.1.2 the power law component 
never seems to disappear completely during the soft States 
A and B. This weak, hard component could be associated 
with the ejection of a small amount of material during State 
A and B, which could explain the low level radio emission 
(also taking into account the possible cooling by soft pho¬ 
tons as explained above). In that case the hardness of the 
spectrum is directly proportional with the strength of the 
radio emission. Although this is also suggested by Meier, 
Koide & Uchida 2001 (see Section 8 . 6 ), we do not find a sig¬ 
nificant relation between the radio flux and the power law 
index, or between the radio flux and the size of the corona 
(estimated by the inner disc radius). 


8.4 Plateau states: the jets and flares 

In Section we discussed the relation between the hard 
State C intervals and the radio oscillation events. We can 
also extend this scenario to describe the presence of the 
quasi-continuous plateau jets, using the argument of contin¬ 
uous ejection during State C. In principle the plateau states 
are just larger, more stable versions of the State C intervals, 
only now a flat spectrum, quasi-continuous jet is observed 
as the jet is being fed from the corona for a longer period 
of time. This allows it to expand further out and hence be¬ 
come optically thin (see also Fig. The combined optically 
thick (at the base) and thin (further out) regions could pro¬ 
duce the flat or sligh^ inverted radio spectrum, as already 
explained in Section & 

In Fig. 1^ we show the radio - infrared spectra for the 
plateau jets and the oscillation events. The similarity in the 
overall shape of the two spectra infers a possible common 
physical origin. 

Two things remain unclear: why does the radio emission 
in class x seems to saturate at a level of ~ 100 mjy, and 
why is the radio emission lower during class X 2 compared 
to the plateau states? As far as t he lo wer radio emission is 
concerned, as explained in Section 5.1 X 2 seems to be differ¬ 


ent from xi and xs in several ways. Although more analysis 
is certainly necessary, the lower X-ray count rate during X 2 
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might reflect an insuflicient supply of material that can be 
injected into the jet. With respect to the saturation level, an 
interaction between the continuous ejection and the possible 
cooling by soft photons might be responsible. 

Summarizing, we have encountered three types of ejec¬ 
tion events (see also Fig. ^): 

• the radio oscillation events, the plateau jets (class x) 
and the low-level radio emission during short State Cs, which 
seem to be related to the same type of mechanism operating 
on different scales; 

• the low-level radio emission during State A and B in¬ 
tervals, which can be related to a hard component but could 
also be a relic from a previous large radio event; 

• the pre- and post-plateau flares. 

The pre- and post-plateau flares differ from the other two 
types of ejection events in several ways. First of all the 
plateau flares are short discrete ejection events which be¬ 
come optically thin very rapidly and, in the case of X 3 sae 
tracked to a distance of ^ 3000 AU from the core (Mirabel 
& Rodriguez 1994, Fender et al. 1999b, Dhawan et al. 2000). 
Secondly, contrary to the low-level radio emission, the radio 
oscillation events and the plateau jets (which are assumed to 
be related to quasi-continuous ejections) the plateau flares 
show clear indications of superluminal motion as separate 
blobs are observed to be ejected. Finally the pre- and post¬ 
plateau flares reach high flux levels up to ~600 mjy (Fender 
et al. 1999b), on average a factor ~3 - 6 higher compared 
to the plateau jets and the radio oscillation events. 

Based on the differences between the plateau flares and 
the other types of radio emissions, we speculate below that 
the plateau flares originate from a different mechanism. 


8.4-1 Plateau flares; a eomparison with state transitions 
in other sources 

In the persistent black hole candidates Cyg X-1 and 
GX 339—4 (Hjellming & Han 1995, Brocksopp et al. 1999, 
Corbel et al. 2000) and in many Soft X-ray Transient Sys¬ 
tems (SXTs) optically thin radio flares have been associated 
with state transitions (Fender & Kuulkers 2001). In the case 
of the SXTs the flares are found during the fast transition 
from the ‘Off’ (quiescence) state to the Soft state, but never 
in the opposite direction (note that the transition from Soft 
to Off in SXTs is usually much slower). For Cyg X-1 and 
possibly also for GX 339—4 the flares have been associated 
with both the Soft - Hard and the Hard - Soft transitions 
(Zhang et al. 1997). The pre-plateau and the post-plateau 
flares, observed in GRS 1915-1-105 could be very closely re¬ 
lated to the flares seen in Cyg X-1, GX 339—4 and the SXTs: 
they occur as the source makes the transition into, or out 
of the plateau state which relates them to Soft - Hard and 
Hard - Soft transitions. 

The relation between the state transitions and the 
optically thin radio flares seems to be a general feature 
for the black hole systems. If the pre- and post-plateau 
flares are also related to state transitions, as argued above, 
GRS 1915-1-105 is one of the few systems to show the flares 
related to both the Soft - Hard and Hard - Soft transi¬ 
tions (note that this indirectly relates these flares to the 
hard State C intervals). In other systems the coverage is of¬ 


ten poor and only occasionally (as in the case of the SXTs) 
flares are observed during the state transitions. 

8.5 Xray : radio correlations during the x-states; 
a comparison to the Low/Hard states in 
GX 339-4 and Cyg X-1 

In two other black hole jet sources, GX 339—4 and Cyg X-1, 
there seems to be a strong three-way correlation between the 
radio, soft and hard X-ray emission during the Low/Hard 
state (Corbel et al. 2000, Brocksopp et al. 1999). During the 
High/Soft states in these two sources, both the radio emis¬ 
sion and the hard X-rays are low. In the Low/Hard state 
GX 339—4 and Cyg X-1 have a radio spectrum which is 
slightly inverted or flat (q~ 0.14 and ~ 0 respectively; see 
Fender 2001 and references therein), and during Low/Hard 
- High/Soft state transitions the sources probably show dis¬ 
crete ejections of expanding relativistic plasma (Corbel et 
al. 2000, Brocksopp et al. 1999). 

With respect to the Xray : radio coupling as presented 
in Section ^ it is clear that during class X 2 GRS 1915-1-105 
resembles both GX 339—4 and Cyg X-1: in all these sources 
there is a correlation between the radio and the soft X-ray 
emission. During xi and X 3 , however, GRS 1915 -fl0 5 shows 
an anticorrelation instead of a correlation (Table js^i]) . As far 
as the radio behaviour is concerned the jets and flares asso¬ 
ciated to xi and xs are most likely also found in GX 339—4 
and Cyg X-1. Although X 2 shows clearly a different radio 
behaviour, with no large radio jets or flares ever observe^ 
the fact remains that the radio spectrum is flat (see Fig 
a flat radio spectrum is considered to be the signature of a 
quasi-continuous jet) as also found for xi, X 35 GX 339—4 
and Cyg X-1. 

Concluding we find that, although there are similarities, 
GRS 1915-1-105 during class x does not completely match 
the behaviour of normal black hole sources in the Low/Hard 
states. 


8.6 Magneto-hydrodynamic jet formation 

The formation of the jets is a process which is probably 
related to magnetohydrodynamic (MHD) effects (see Meier 
et al. (2001) for a review). In that case, the jets are pro¬ 
duced/confined by the toroidal coiling of the magnetic held 
lines (which are embedded in the accretion disc) caused by 
the differential rotation of the accreting material (Meier et 
al. 2001). As stronger magnetic fields are expected for ge¬ 
ometrically thick accretion flows, the appearance of a jet 
is directly linked to the Low/Hard states of black hole and 
neutron star X-ray binary systems (geometrically thick flows 
such as an ADAF or a corona are believed to be responsi¬ 
ble for the production of the hard component in the spec¬ 
trum). In fact, Meier (2001) shows that for X-ray binaries 
the most powerful jets are only produced in geometrically 
thick inner accretion flows where strong vertical (poloidal) 
magnetic fields are induced which are able to produce pow¬ 
erful (~10®® erg s“^) radio jets. During the High/Soft states, 
on the other hand, when the spectrum is dominated by the 
soft multi-temperature black body component of the geo¬ 
metrically thin disc, the power in the jet is predicted to be 
a hundred times weaker compared to the thick flow solution 
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(Meier 2001). This is in agreement with both our results 
and the observed relation between the black hole state and 
the jet in GX 339—4 by Fender et al. (1999a), and it could 
explain why during the soft States A and B only low-level 
radio emission is observed. 

Although other processes such as shock fronts forming 
in the jets have been suggested and MHD processes still 
have difficulty explaining some of the high outflow speeds 
in the jet systems (see Meier et al. 2001) we feel, based on 
the arguments given above, that MHD processes are a good 
candidate for explaining the formation of the jets in the X- 
ray binary systems. Note however, that Meier (2001) and 
Meier et al. (2001) relate the strength of the radio emission 
to the hardness o f th e X-ray spectrum, while in our simple 
scenario (Section B.l) it is related to the length of the State 
C intervals. 


8.7 Alternatives 

In papers by Naik et al. (1999), Naik & Rao (2000), Eiken- 
berry et al. (2000) and Yadav (2001) it is concluded that that 
soft dips are related to the radio oscillation events. The ob¬ 
servations presented in the papers can be classified as class 
/3, A, 6 or V. We therefore suggest that in these observa¬ 
tions the long (^ 100s), hard State C dips, which are always 
present just prior to State A dips (in the case when both 
States A and C are present; Belloni et al. 2000a), are re¬ 
sponsible for the production of the radio oscillation events 
as explained in the previous sections. The fact that we can 
come to this conclusion is related to our much larger X-ray 
- radio overlap compared to that of the authors mentioned 
above. 

Naik et al. (1999) and Naik & Rao (2000) also suggest 
that the combination of several soft dips produces the large 
pre- and post-plateau flares. Although we recognize that sev¬ 
eral small ejections can be combined to give one big radio 
event (such as a flare), we do not believe that is the case 
here. Probably the most convincing argument against this 
is the fact that we see several observations with repeating 
hard X-ray dips (e^Figs. |3a,^ and ^) or with repeating 
soft dips (eg. Fig. and m were soft State A dips are 
preceded by State C dips), during which no large plateau 
flares are reported (note that the argument presented here 
does not depend on whether the dips are hard or soft: in the 
model of Naik et al. (1999) and Naik & Rao (2000) the soft 
dips are responsible for the ejection of material, in our model 
the hard dips are responsible). Moreover, Pooley & Fender 
(1997) and Fender & Pooley (1998, 2000) report multiple oc¬ 
casions of long sequences of radio oscillation events with no 
major flare. Finally, a close examination of similar outbursts 
to the one discussed by Naik et al. (1999) reveals that the 
hrst series of long State C dips (observations K-52-01 and 
K-52-02, see Table |^), and therefore also the corresponding 
radio oscillation events, start after the large jet is observed 
(see figure 4 in Fender et al. 1999b). Taking into account the 
delay between the X-rays and the radio emission, it is very 
unlikely that the combined effect of these small ejections 
could be responsible for the production of the large plateau 
jets. 

In two different papers Nayakshin, Rappaport, Melia 
(2000) and Janiuk, Czerny, Siemiginowska (2000) conclude 
they need an outflow in the form of a jet, to correctly model 


Table 4. The emipircal X-ray : radio correlation in 
GRS 1915-1-105. 


Radio Behaviour 

X-ray behaviour 

Radio oscillation events 

repeating, hard, well 
separated ^ 100s State C 

Low level emission 

adjacent, ^ 100s 

State C or A, B 

Plateau Jet 

long uninterrupted 


State C 

Pre-, Post-plateau flares 

State transitions? 


the complex behaviour of GRS 1915-1-105. They both model 
the jet outflow in the same way assuming the ejections set in 
as the source makes the transition into the High/Soft state. 
As explained in the previous sections, this is inconsistent 
with our results: we suggested, based observational facts, 
that the ejection of material does not occur at a single point 
but is more or less continuous during State C. Therefore 
we feel that future models should incorporate jet outflows 
which are related to hard (State C) intervals which could, 
in principle, lead to the production of the radio oscillation 
events and the plateau jet in a natural way (as explained in 
Section 8.1). 


9 SUMMARY & CONCLUSIONS 

This work has clearly established that radio emission from 
GRS 1915-1-105 is intimately related to the presence of hard 
(power-law-dominated) intervals in the X-ray light curves. 
This in turn physically implies a clear relation between a 
radiatively inefficient flow close to the black hole, and a 
synchrotron-emitting outflow or jet. 

In Table ^ we summerize the empirical X-ray : radio 
correlation, moreover we hnd that: 

• The (high flux) radio oscillation events are closely re¬ 
lated to a series of long (or low count rate), well-separated, 
spectrally hard State C intervals in the X-ray light curve. 

• The State G intervals must be at least 100s in length 
and have a recurrence time that allows for the individual 
radio oscillation events to be observed. 

• We conhrm a one-to-one relation between radio oscilla¬ 
tion events and series of long State C intervals: each State 
C interval produces a new “radio flare”. 

• Short (Si 100s) State G intervals and State A and B 
intervals produce a similar low-level radio emission. 

• The radio oscillation profiles seem to be determined by 
particle injection/acceleration (rise phase) and adiabatic ex¬ 
pansion losses (decay phase), which are both frequency in¬ 
dependent effects. 

• During class xi and xs (plateau states) we find a quasi- 
continuous jet with a flat synchrotron spectrum (extending 
to at least the near-infrared) and a flux level of ^ 100 mjy, 
an anticorrelation between the radio flux and the soft (ASM) 
X-ray flux, and a correlation between the radio flux and the 
ASM HRa colour. 

• During class xz we also find a flat synchrotron spec¬ 
trum, but a flux level which is about a factor 10 lower, a 












Hard X-ray states and radio emission in GRS 1915-1-105 21 


correlation between the radio flux and the soft (ASM) X- 
ray flux and a anticorrelation between the radio flux and 
the ASM HR 2 colour. 

Based on this we suggest a simple scenario which 
describes the complex X-ray : radio behaviour of 
GRS 1915-tl05: 

• During State C a more or less continuous ejection of 
relativistic particles takes place; 

• The length (or depth, as longer also means deeper) of 
the State C interval determines the strength of the radio 
emission: longer State Cs have a higher radio flux; 

• The separation between the State Cs determines the 
shape of the radio light curve: if the intervals are to close 
together one observes a “stable” radio flux level, while if 
the separation is large enough for the individual radio flares 
to be observed one finds a radio light curve dominated by 
oscillations. 

Within the simple scenario we can explain: 

• the radio oscillation events: series of well separated, long 
State Cs; 

• the radio plateau jets: uninterrupted State Cs; 

• the low-level radio emission during the short State Cs: 
short, near-contiguous State Cs. 

All these types of radio emission seem to be related to the 
same physical mechanism operating on different scales. 

However, this scenario does not explain: 

• The disproportional higher radio flux in the case of the 
long State C compared to the short State C intervals (al¬ 
though cooling by soft X-ray photons might be responsible); 

• The low-level radio emission in the case of State A and 
B, which cannot be explained based on argument of a con¬ 
tinuous ejection during State C; 

• The radio plateau flares observed prior to and just af¬ 
ter the plateau states, which seem to be related to the Soft 
- Hard and Hard - Soft X-ray state transitions as also ob¬ 
served in SXTs; 

• The lower radio emission during class X 2 compared to 
classes xi and xs- 

With respect to the low-level radio emission found dur¬ 
ing States A, B and short State C intervals, we cannot rule 
out the possibility that the emission is a relic from a previous 
large radio event. 

The most important result is probably the establish¬ 
ment of the State C - radio relation, and in particu¬ 
lar the suggested one-to-one relation between the long 
State C intervals and the radio oscillation events. Because 
GRS 1915-1-105 has a much larger coverage in the radio than 
in the X-rays (PCA observations), this intimate relation can 
provide us with much more information about the behaviour 
of the source its behaviour than is revealed to us by the X- 
rays alone. However, due to the time delay between the radio 
and the X-rays, the disc - jet connection is probably better 
studied in the X-ray and IR bands. 

The analysis presented here, together with previous re¬ 
ports on other black hole X-ray binary systems such as 
GX 339—4 (Fender et al. 1999a) and XTE J1550-564 (Cor¬ 
bel et al. 2001 ), provides strong evidence that hard, power 
law-dominated states can support a radio-emitting outflow. 


This suggest, as discussed in Section S.6, that MHD effects 
could be responsible for the production and/or confinement 
of the jets found in these systems. 
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